Improving service efficiency and throughput of cardiac surgery patients using Monte Carlo simulation: a queueing setting.

Bed occupancy rate (BOR) is important for healthcare policymakers. Studies showed the necessity of using simulation approach when encountering complex real-world problems to plan the optimal use of resources and improve the quality of services. So, the aim of the present study is to estimate average length of stay (LOS), BOR, bed blocking probability (BBP), and throughput of patients in a cardiac surgery department (CSD) using simulation models. We studied the behavior of a CSD as a complex queueing system at the Farshchian Hospital. In the queueing model, customers were patients and servers were beds in intensive care unit (ICU) and post-operative ward (POW). A computer program based on the Monte Carlo simulation, using Python software, was developed to evaluate the behavior of the system under different number of beds in ICU and POW. The queueing simulation study showed that, for a fixed number of beds in ICU, BOR in POW decreases as the number of beds in POW increases and LOS in ICU increases as the number of beds in POW decreases. Also, based on the available data, the throughput of patients in the CSD during 800 days was 1999 patients. Whereas, the simulation results showed that, 2839 patients can be operated in the same period. The results of the simulation study clearly demonstrated the behavior of the CSD; so, it must be mentioned, hospital administrators should design an efficient plan to increase BOR and throughput of patients in the future.

Cortical thickness analysis in temporal lobe epilepsy using fully Bayesian spectral method in magnetic resonance imaging.

BACKGROUND: Temporal lobe epilepsy (TLE) is the most common type of epilepsy associated with changes in the cerebral cortex throughout the brain. Magnetic resonance imaging (MRI) is widely used for detecting such anomalies; nevertheless, it produces spatially correlated data that cannot be considered by the usual statistical models. This study aimed to compare cortical thicknesses between patients with TLE and healthy controls by considering the spatial dependencies across different regions of the cerebral cortex in MRI. METHODS: In this study, T1-weighted MRI was performed on 20 healthy controls and 33 TLE patients. Nineteen patients had a left TLE and 14 had a right TLE. Cortical thickness was measured for all individuals in 68 regions of the cerebral cortex based on images. Fully Bayesian spectral method was utilized to compare the cortical thickness of different brain regions between groups. Neural networks model was used to classify the patients using the identified regions. RESULTS: For the left TLE patients, cortical thinning was observed in bilateral caudal anterior cingulate, lateral orbitofrontal (ipsilateral), the bilateral rostral anterior cingulate, frontal pole and temporal pole (ipsilateral), caudal middle frontal and rostral middle frontal (contralateral side). For the right TLE patients, cortical thinning was only observed in the entorhinal area (ipsilateral). The AUCs of the neural networks for classification of left and right TLE patients versus healthy controls were 0.939 and 1.000, respectively. CONCLUSION: Alteration of cortical gray matter thickness was evidenced as common effect of epileptogenicity, as manifested by the patients in this study using the fully Bayesian spectral method by taking into account the complex structure of the data.

Factors affecting systolic blood pressure trajectory in low and high activity conditions.

Background: Typically, blood pressure dips during sleep and increases during daytime. The blood pressure trend is affected by the autonomic nervous system. The activity of this system is observable in the low and high activity conditions. The aim of this study was to assess the effect of individual characteristics on systolic blood pressure (SBP) across day-night under low and high activity conditions. Methods: The samples were 34 outpatients who were candidates for evaluation of 24 hours of blood pressure with an ambulatory. They were admitted to the heart clinic of Farshchian hospital, located in Hamadan province in the west of Iran. The hourly SBP during 24 hours was considered as a response variable. To determine the factors effecting SBP in each condition, the hidden semi-Markov model (HSMM), with 2 hidden states of low and high activity, was fitted to the data. Results: Males had lower SBP than females in both states. The effect of age was positive in the low activity state (ß=0.30; p<0.001) and negative in high activity state (ß= -0.21; p=0.001). The positive effect of cigarette smoking on SBP was seen in low activity state (ß=5.02; p=0.029). The overweight and obese patients had higher SBP compared to others in high activity state (ß=11.60; p<0.001 and ß=5.87; p=0.032, respectively). Conclusion: The SBP variability can be displayed by hidden states of low and high activity. Moreover, the effects of studied variables on SBP were different in low and high activity states.

Regularized Weighted Nonparametric Likelihood Approach for High-Dimension Sparse Subdistribution Hazards Model for Competing Risk Data.

Variable selection and penalized regression models in high-dimension settings have become an increasingly important topic in many disciplines. For instance, omics data are generated in biomedical researches that may be associated with survival of patients and suggest insights into disease dynamics to identify patients with worse prognosis and to improve the therapy. Analysis of high-dimensional time-to-event data in the presence of competing risks requires special modeling techniques. So far, some attempts have been made to variable selection in low- and high-dimension competing risk setting using partial likelihood-based procedures. In this paper, a weighted likelihood-based penalized approach is extended for direct variable selection under the subdistribution hazards model for high-dimensional competing risk data. The proposed method which considers a larger class of semiparametric regression models for the subdistribution allows for taking into account time-varying effects and is of particular importance, because the proportional hazards assumption may not be valid in general, especially in the high-dimension setting. Also, this model relaxes from the constraint of the ability to simultaneously model multiple cumulative incidence functions using the Fine and Gray approach. The performance/effectiveness of several penalties including minimax concave penalty (MCP); adaptive LASSO and smoothly clipped absolute deviation (SCAD) as well as their L2 counterparts were investigated through simulation studies in terms of sensitivity/specificity. The results revealed that sensitivity of all penalties were comparable, but the MCP and MCP-L2 penalties outperformed the other methods in term of selecting less noninformative variables. The practical use of the model was investigated through the analysis of genomic competing risk data obtained from patients with bladder cancer and six genes of CDC20, NCF2, SMARCAD1, RTN4, ETFDH, and SON were identified using all the methods and were significantly correlated with the subdistribution.

Predicting the incidence of brucellosis in Western Iran using Markov switching model.

OBJECTIVE: Brucellosis is a zoonosis almost chronic disease. Brucellosis bacteria can remain in the environment for a long time. Thus, climate irregularities could pave the way for the survival of the bacterium brucellosis. Brucellosis is more common in men 25 to 29 years of age, in the western provinces, and in the spring months. The aim of this study is to investigate the effect of climatic factors as well as predicting the incidence of brucellosis in Qazvin province using the Markov switching model (MSM). This study is a secondary study of data collected from 2010 to 2019 in Qazvin province. The data include brucellosis cases and climatic parameters. Two state MSM with time lags of 0, 1 and 2 was fitted to the data. The Bayesian information criterion (BIC) was used to evaluate the models. RESULTS: According to the BIC, the two-state MSM with a 1-month lag is a suitable model. The month, the average-wind-speed, the minimum-temperature have a positive effect on the number of brucellosis, the age and rainfall have a negative effect. The results show that the probability of an outbreak for the third month of 2019 is 0.30%.

Generalized linear mixed hidden semi-Markov models in longitudinal settings: A Bayesian approach.

Hidden Markov and semi-Markov models (H(S)MMs) constitute useful tools for modeling observations subject to certain dependency structures. The hidden states render these models very flexible and allow them to capture many different types of latent patterns and dynamics present in the data. This has led to the increased popularity of these models, which have been applied to a variety of problems in various domains and settings, including longitudinal data. In many longitudinal studies, the response variable is categorical or count-type. Generalized linear mixed models (GLMMs) can be used to analyze a wide range of variables, including categorical and count. The present study proposes a model that combines HSMMs with GLMMs, leading to generalized linear mixed hidden semi-Markov models (GLM-HSMMs). These models can account for time-varying unobserved heterogeneity and handle different response types. Parameter estimation is achieved using a Monte Carlo Newton-Raphson (MCNR)-like algorithm. In our proposed model, the distribution of the random effects depends on hidden states. We illustrate the applicability of GLM-HSMMs with an example in the field of occupational health, where the response variable consists of count values. Furthermore, we assess the performance of our MCNR-like algorithm through a simulation study.

Predicting preeclampsia and related risk factors using data mining approaches: A cross-sectional study.

BACKGROUND: Preeclampsia is a type of pregnancy hypertension disorder that has adverse effects on both the mother and the fetus. Despite recent advances in the etiology of preeclampsia, no adequate clinical screening tests have been identified to diagnose the disorder. OBJECTIVE: We aimed to provide a model based on data mining approaches that can be used as a screening tool to identify patients with this syndrome and also to identify the risk factors associated with it. MATERIALS AND METHODS: The data used to perform this cross-sectional study were extracted from the clinical records of 726 mothers with preeclampsia and 726 mothers without preeclampsia who were referred to Fatemieh Hospital in Hamadan City during April 2005-March 2015. In this study, six data mining methods were adopted, including logistic regression, k-nearest neighborhood, C5.0 decision tree, discriminant analysis, random forest, and support vector machine, and their performance was compared using the criteria of accuracy, sensitivity, and specificity. RESULTS: Underlying condition, age, pregnancy season and the number of pregnancies were the most important risk factors for diagnosing preeclampsia. The accuracy of the models were as follows: logistic regression (0.713), k-nearest neighborhood (0.742), C5.0 decision tree (0.788), discriminant analysis (0.687), random forest (0.758) and support vector machine (0.791). CONCLUSION: Among the data mining methods employed in this study, support vector machine was the most accurate in predicting preeclampsia. Therefore, this model can be considered as a screening tool to diagnose this disorder.

The Impact of Control Measures and Holiday Seasons on Incidence and Mortality Rate of COVID-19 in Iran.

BACKGROUND: Preventive measures on the COVID-19 pandemic is an effective way to control its spread. We aimed to investigate the effect of control measures and holiday seasons on the incidence and mortality rate of COVID-19 in Iran. STUDY DESIGN: An observational study. METHODS: The daily data of confirmed new cases and deaths in Iran were taken from the Johns Hopkins University COVID-19 database. We calculated weekly data from 19 Feb to 6 Oct 2020. To estimate the impact of control measures and holiday seasons on the incidence rate of new cases and deaths, an autoregressive hidden Markov model (ARHMM) with two hidden states fitted the data. The hidden states of the fitted model can distinguish the peak period from the non-peak period. RESULTS: The control measures with a delay of one-week and two-week had a decreasing effect on the new cases in the peak and non-peak periods, respectively (P=0.005). The holiday season with a two-week delay increased the total number of new cases in the peak periods (P=0.031). The peak period for the occurrence of COVID-19 was estimated at 3 weeks. In the peak period of mortality, the control measures with a three-week delay decreased the COVID-19 mortality (P=0.010). The expected duration of staying in the peak period of mortality was around 6 weeks. CONCLUSION: When an increasing trend was seen in the country, the control measures could decline the incidence and mortality related to COVID-19. Implementation of official restrictions on holiday seasons could prevent an upward trend of incidence for COVID-19 during the peak period.

Multistate recursively imputed survival trees for time-to-event data analysis: an application to AIDS and mortality post-HIV infection data.

BACKGROUND: This study aimed to introduce recursively imputed survival trees into multistate survival models (MSRIST) to analyze these types of data and to identify the prognostic factors influencing the disease progression in patients with intermediate events. The proposed method is fully nonparametric and can be used for estimating transition probabilities. METHODS: A general algorithm was provided for analyzing multi-state data with a focus on the illness-death and progressive multi-state models. The model considered both beyond Markov and Non-Markov settings. We also proposed a multi-state random survival method (MSRSF) and compared their performance with the classical multi-state Cox model. We applied the proposed method to a dataset related to HIV/AIDS patients based on a retrospective cohort study extracted in Tehran from April 2004 to March 2014 consist of 2473 HIV-infected patients. RESULTS: The results showed that MSRIST outperformed the classical multistate method using Cox Model and MSRSF in terms of integrated Brier score and concordance index over 500 repetitions. We also identified a set of important risk factors as well as their interactions on different states of HIV and AIDS progression. CONCLUSIONS: There are different strategies for modelling the intermediate event. We adapted two newly developed data mining technique (RSF and RIST) for multistate models (MSRSF and MSRIST) to identify important risk factors in different stages of the diseases. The methods can capture any complex relationship between variables and can be used as a useful tool for identifying important risk factors in different states of this disease.

Relationship between Urbanization and Cancer Incidence in Iran Using Quantile Regression.

Quantile regression is an efficient method for predicting and estimating the relationship between explanatory variables and percentile points of the response distribution, particularly for extreme percentiles of the distribution. To study the relationship between urbanization and cancer morbidity, we here applied quantile regression. This cross-sectional study was conducted for 9 cancers in 345 cities in 2007 in Iran. Data were obtained from the Ministry of Health and Medical Education and the relationship between urbanization and cancer morbidity was investigated using quantile regression and least square regression. Fitting models were compared using AIC criteria. R (3.0.1) software and the Quantreg package were used for statistical analysis. With the quantile regression model all percentiles for breast, colorectal, prostate, lung and pancreas cancers demonstrated increasing incidence rate with urbanization. The maximum increase for breast cancer was in the 90th percentile (ß=0.13, p-value<0.001), for colorectal cancer was in the 75th percentile (ß=0.048, p-value<0.001), for prostate cancer the 95th percentile (ß=0.55, p-value<0.001), for lung cancer was in 95th percentile (ß=0.52, p-value=0.006), for pancreas cancer was in 10th percentile (ß=0.011, p-value<0.001). For gastric, esophageal and skin cancers, with increasing urbanization, the incidence rate was decreased. The maximum decrease for gastric cancer was in the 90th percentile(ß=0.003, p-value<0.001), for esophageal cancer the 95th (ß=0.04, p-value=0.4) and for skin cancer also the 95th (ß=0.145, p-value=0.071). The AIC showed that for upper percentiles, the fitting of quantile regression was better than least square regression. According to the results of this study, the significant impact of urbanization on cancer morbidity requirs more effort and planning by policymakers and administrators in order to reduce risk factors such as pollution in urban areas and ensure proper nutrition recommendations are made.

Predicting the Survival Time for Bladder Cancer Using an Additive Hazards Model in Microarray Data.

BACKGROUND: One substantial part of microarray studies is to predict patients' survival based on their gene expression profile. Variable selection techniques are powerful tools to handle high dimensionality in analysis of microarray data. However, these techniques have not been investigated in competing risks setting. This study aimed to investigate the performance of four sparse variable selection methods in estimating the survival time. METHODS: The data included 1381 gene expression measurements and clinical information from 301 patients with bladder cancer operated in the years 1987 to 2000 in hospitals in Denmark, Sweden, Spain, France, and England. Four methods of the least absolute shrinkage and selection operator, smoothly clipped absolute deviation, the smooth integration of counting and absolute deviation and elastic net were utilized for simultaneous variable selection and estimation under an additive hazards model. The criteria of area under ROC curve, Brier score and c-index were used to compare the methods. RESULTS: The median follow-up time for all patients was 47 months. The elastic net approach was indicated to outperform other methods. The elastic net had the lowest integrated Brier score (0.137±0.07) and the greatest median of the over-time AUC and C-index (0.803±0.06 and 0.779±0.13, respectively). Five out of 19 selected genes by the elastic net were significant (P<0.05) under an additive hazards model. It was indicated that the expression of RTN4, SON, IGF1R and CDC20 decrease the survival time, while the expression of SMARCAD1 increase it. CONCLUSION: The elastic net had higher capability than the other methods for the prediction of survival time in patients with bladder cancer in the presence of competing risks base on additive hazards model.

Competing risks data analysis with high-dimensional covariates: an application in bladder cancer.

Analysis of microarray data is associated with the methodological problems of high dimension and small sample size. Various methods have been used for variable selection in high-dimension and small sample size cases with a single survival endpoint. However, little effort has been directed toward addressing competing risks where there is more than one failure risks. This study compared three typical variable selection techniques including Lasso, elastic net, and likelihood-based boosting for high-dimensional time-to-event data with competing risks. The performance of these methods was evaluated via a simulation study by analyzing a real dataset related to bladder cancer patients using time-dependent receiver operator characteristic (ROC) curve and bootstrap .632+ prediction error curves. The elastic net penalization method was shown to outperform Lasso and boosting. Based on the elastic net, 33 genes out of 1381 genes related to bladder cancer were selected. By fitting to the Fine and Gray model, eight genes were highly significant (P<0.001). Among them, expression of RTN4, SON, IGF1R, SNRPE, PTGR1, PLEK, and ETFDH was associated with a decrease in survival time, whereas SMARCAD1 expression was associated with an increase in survival time. This study indicates that the elastic net has a higher capacity than the Lasso and boosting for the prediction of survival time in bladder cancer patients. Moreover, genes selected by all methods improved the predictive power of the model based on only clinical variables, indicating the value of information contained in the microarray features.

Predicting hepatitis B monthly incidence rates using weighted Markov chains and time series methods.

BACKGROUND: Hepatitis B (HB) is a major global mortality. Accurately predicting the trend of the disease can provide an appropriate view to make health policy disease prevention. This paper aimed to apply three different to predict monthly incidence rates of HB. METHODS: This historical cohort study was conducted on the HB incidence data of Hamadan Province, the west of Iran, from 2004 to 2012. Weighted Markov Chain (WMC) method based on Markov chain theory and two time series models including Holt Exponential Smoothing (HES) and SARIMA were applied on the data. The results of different applied methods were compared to correct percentages of predicted incidence rates. RESULTS: The monthly incidence rates were clustered into two clusters as state of Markov chain. The correct predicted percentage of the first and second clusters for WMC, HES and SARIMA methods was (100, 0), (84, 67) and (79, 47) respectively. CONCLUSIONS: The overall incidence rate of HBV is estimated to decrease over time. The comparison of results of the three models indicated that in respect to existing seasonality trend and non-stationarity, the HES had the most accurate prediction of the incidence rates.

High-dimensional additive hazards regression for oral squamous cell carcinoma using microarray data: a comparative study.

Microarray technology results in high-dimensional and low-sample size data sets. Therefore, fitting sparse models is substantial because only a small number of influential genes can reliably be identified. A number of variable selection approaches have been proposed for high-dimensional time-to-event data based on Cox proportional hazards where censoring is present. The present study applied three sparse variable selection techniques of Lasso, smoothly clipped absolute deviation and the smooth integration of counting, and absolute deviation for gene expression survival time data using the additive risk model which is adopted when the absolute effects of multiple predictors on the hazard function are of interest. The performances of used techniques were evaluated by time dependent ROC curve and bootstrap .632+ prediction error curves. The selected genes by all methods were highly significant (P < 0.001). The Lasso showed maximum median of area under ROC curve over time (0.95) and smoothly clipped absolute deviation showed the lowest prediction error (0.105). It was observed that the selected genes by all methods improved the prediction of purely clinical model indicating the valuable information containing in the microarray features. So it was concluded that used approaches can satisfactorily predict survival based on selected gene expression measurements.